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ABSTRACT 

Void regions of the Universe offer a special environment for studying cosmology 
and galaxy formation, which may expose weaknesses in our understanding of these 
phenomena. Although galaxies in voids are observed to be predominately gas rich, 
star forming and blue, a sub-population of bright red void galaxies can also be found, 
whose star formation was shut down long ago. Are the same processes that quench 
star formation in denser regions of the Universe also at work in voids? 

We compare the luminosity function of void galaxies in the 2dF Galaxy Redshift 
Survey, to those from a galaxy formation model built on the Millennium Simulation. 
We show that a global star formation suppression mechanism in the form of low lumi- 
nosity "radio mode" AGN heating is sufhcient to reproduce the observed population 
of void early-types. Radio mode heating is environment independent other than its 
dependence on dark matter halo mass, where, above a critical mass threshold of ap- 
proximately Mvir ~ 10^^-^ Mq, gas cooling onto the galaxy is suppressed and star 
formation subsequently fades. In the Millennium Simulation, the void halo mass func- 
tion is shifted with respect to denser environments, but still maintains a high mass 
tail above this critical threshold. In such void halos, radio mode heating remains effi- 
cient and red galaxies are found; collectively these galaxies match the observed space 
density without any modification to the model. Consequently, galaxies living in vastly 
different large-scale environments but hosted by halos of similar mass are predicted to 
have similar properties, consistent with observations. 

Key words: galaxies: statistics, galaxies: evolution, galaxies: active, large-scale struc- 
ture of the Universe 



1 INTRODUCTION 

According to conventional cosmology, voids provide a 
special laboratory for studying galaxy formation, in which 
the dark matter density is locally that of a low den- 
sity Universe. Peebles has claimed l|Peeblesl 1200 ll ) that 
A cold dark matter (ACDM) simulations overproduce 
void galaxies in comparison to observation and he sug- 
gested that this may be evidence of as-yet-undiscovered 
fundamental physics not included in ACDM cosmology. 
Previous work has laid the grou ndwork for examining 



this question quantitatively (e.E. iMathis & White 




2002; 


Benson et al.1 2003bl: Mo et al. 2004: Croton et al 




2004, 


20051: iHovle & Vogeley"2004"; 'Hovlc ot al."2005'; 'Roias 


et al. 


2004. [2005: Colbcrg et al. 2005; Goldberg et al] 




20051: 


Clemens et al.1 20061: IXia et al. bood: IPatiri et al. 




2009: 


Sorrentino et al.l l2006l: Ivon Benda-Beckmann & MiilleiJ 


20081). In Croton et al.l (|2005f) we introduced a oractical 



definition of the void luminosity function and measured it 



for early- and late-type galaxies in the 2dF Galaxy Redshift 
Survey (2dFGRS). Here, we investigate whether ACDM 
simulations can describe the observed void luminosity 
functions quantitatively, or whether new physics may be 
required. 

State-of-the-art modelling is not sufficiently advanced 
to treat all the physical processes of galaxy formation from 
first principles. Here we adopt a semi-analytic approach 
(|White fc Frenklll99ll) in which the Millennium Simulation 
of dark matter is coupled with an analytic simulation for 
baryonic evolution in a cosmologically evolving background. 
A simplifying approximation made in the current genera- 
tion of semi-analytic models is that the process of galaxy 
formation depends only indirectly on the environment or 
age of a given dark matter halo, via their impact on the 
halo mass function. Wh ile this clearly cannot be exactly true 
(e.g. lCroton et aDl2007h . it is interesting to ask whether such 
simplified models are capable of reproducing the observed 
luminosity distributions of early- and late-type galaxies in 
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voids within the observational uncertainties. As we show in 
the following, the answer is affirmative. 

Key to the success of semi-analytic models in describ- 
ing the global galaxy luminosity function is t he presence of 
heatin g mechanisms to shut off star formation l|Benson et al.l 
l2003al ). The exact processes operating are currently subject 
to some debate. What is not in dispute is the observational 
evidence of galaxies at the bright-end of the galaxy lumi- 
nosity function with ongoing, low-luminosity active nuclei. 
It is thus no surprise that most current semi-analytic mod- 
els incorpo rate active galactic n uclei (AGN) heating as a key 
ingredient. ICroton et al.l (|2006l ) called this the "radio mode" 
of AGN evolution. Another process which can shut off star 
formation in subhalos is the stripping of gas when a satellite 
falls into a larger halo, called "strangulation". 

There is clear evidence of a population of "red and 
dead" early-type galaxie^ even in voids (ICroton et al.l 
I2OO5I ). yet it is not obvious whether either or both of these 
mechanisms (i.e. AGN heating and strangulation) can oper- 
ate efficiently in voids or whether some other mechanism is 
needed. We investigate this question via semi-analytic sim- 
ulations, and find that strangulation plays a minor role. We 
find that heating when the central galaxy passes the critical 
mass is sufficient to account for the observed abundance of 
early-type galaxies in voids. The different luminosity func- 
tions of early-type galaxies between voids and mean density 
environments can thus be entirely attributed to the differ- 
ence in the halo mass function in the two environments, 
within the accuracy of present measurements. This expla- 
nation for the early-type luminosity function has observ- 
able consequences as discussed below. If these are verified, 
it will lend credence to the physical interpretation of the 
semi-analytic model. 

This paper is organised as follows. Sections [2] and 
[3] describe the galaxy formation model and our measure 
of environment within the Millennium Simulation, respec- 
tively. We present our void galaxy analysis in Section |4] 
and discuss our results in light of recent work and also 
within the broader context of galaxy formation theory. Fi- 
nally, Section [5] provides a brief summary. Throughout we 
assume a standard WMAP first yea r ACDM cosmology 
jSpergel et al.ll2003l : [Seliak et al.ll2005^ and Hubble param- 
eter Ho = 100/i"^kms"^Mpc"^ 



2 THE GALAXY FORMATION MODEL 

The galaxy formation model we use to study void e n viron - 
ments is identical to that described in ICroton et all (|2006l ) 
(including parameter choices). This model is implemented 
on top of the Millennium R un ACDM dark matter simu- 
lation (|Springel et al.l [20051 ). Below we briefiy outline the 
relevant aspects of the simulation and model to our current 
work, and refer the interested reader to the above references 
for further information. 

The Millennium Run simulation follows the dynami- 
cal evolution of 10^^" dark matter particles in a periodic 
box of side- length 500/!.~^Mpc with a mass resolution per 

^ The term "red and dead" refers to those galaxies on the red 
sequence which have had no active star formation in the last 
several Gyr 



particle of 8.6 x 10* /i ^M0. It adopts cosmological pa- 
rameter v alues consistent wit h a combined analysis of the 
2dFGRS (CoUess et al' '200l') and first year WMAP data 
(Sp ergel et al.. ,2003, : .Seliak et al.. ■2005. ') . Friends-of-friends 
(FOF) hales are identified in the simulation using a linking 
length of 0.2 the mean particle separation, while substruc- 
ture within each FOF halo is found with an im proved and ex- 
tcndcd version of the SUBFIND algorithm of ISpringel et al.] 
(2001). Having determined all halos and subhalos at all 
output snapshots, the hierarchical merging trees are con- 
structed; these describe in detail how structures grow as the 
universe evolves. These trees form the backbone onto which 
we couple our model of galaxy formation. 

Inside each tree, virialised dark matter halos at each 
redshift are assumed to attract ambient gas from the sur- 
rounding medium, from which galaxies form and evolve. Our 
model effectively tracks a wide range of galaxy formation 
physics in each halo using simple parametrised forms, in- 
cluding reionization of the inter-galactic medium at high 
redshift, radiative cooling of hot gas and the formation of 
cooling fiows, star formation in the cold disk and the re- 
sulting supernova feedback, black hole growth and AGN 
feedback through the 'quasar' and 'radio' epochs of AGN 
evolution, metal enrichment of the inter-galactic and intra- 
cluster medium, and galaxy morphology shaped through 
mergers and merger-induced starbursts. At z = the galaxy 
formation model contains approximately 9 million galaxies 
brighter than our completeness limit of Mbj — 51ogjQft = 
— 15.8 (representing the mean luminosity of a central galaxy 
in a halo containing 64 particles). 

A critical feature for the success of this model is the 
inclusion of an AGN component and the separation of AGN 
into their high and low accretion states, called the "quasar 
mode" and "radio mode" respectively. Within this picture 
the two AGN modes are di stinct in their cause and effect. 
Importantly for the work of lCroton et al. 

1 diood), the radio 

mode was used to suppress the cooling of gas onto central 
galaxies living in group and cluster-sized halos, and is pri- 
marily important at late times (typically z < 1) and not 
earlier (where the SFR density of the universe peaks, e.g. 
iMadau et al.|[l99^ ). This model was tuned to provide a good 
match to the local luminosity function and colours of galax- 
ies. The assumed radio mode black hole accretion rate, 

mBH oc rriBH /hot Kdr , (1) 

is very efficient above a critical halo mass threshold which is 
approximately constant with time. This mass turns out to 
be Mvir ~ 1O^'^~^''M0. From this point-of-view, halos grow 
with time until they cross the critical halo mass, after which 
the radio mode dominates, cooling gas slows, then stops, and 
star formation shut-down follows. Note that no environment 
dependence has been assumed in this implementation of ra- 
dio mode heating. We will return to this point in Section 



3 ESTIMATING LOCAL DENSITIES 

To facilitate a fair comparison between the model and 
observed luminosity functions we determine environments 
in the Millennium Simulation box in the same way as 
ICroton et all l|2005l ') did for the 2dFGRS. This is actually 
much simpler to do than with a magnitude limited redshift 
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Figure 1. T he galaxy luminosity functions for all environments 
(filled circles. rNorberg _et_ajj200^ and in void regions exclusively 
(open squares. ICroton et al.ll2005l) . both from the 2dFGRS. The 
matching solid lines indicate the semi-analytic model result anal- 
ysed in the same way. The model is tuned for good agreement 
•with the global 2dFGRS luminosity function; from this alone the 
distribution of galaxy magnitudes in voids is also reproduced. 



Figure 2. The galaxy luminosity function in voids, now bro- 
ken into their blue/late-type and red/early- type sub-populations 
(see figure legend). The dotted lines show the model luminosity 
functions for central galaxies only, i.e. when satellites are removed 
(see Section [4.2l l. Both red and blue sequence model galaxies have 
similar luminosity distributions to the 2dFGRS data. The model 
early-type void population is dominated by central galaxies. 



survey, as one has a perfect measure of survey geometry 
(a box in our case with volume Vbox) and selection func- 
tion (uniform and complete). We consider the redshift zero 
snapshot of the model, which we call the "local" population. 
Note that, although the 2dFGRS has a median redshift of 
Z'^0.1, all 2dFGRS magnitudes are k-corrected to z = 0, and 
for our purposes we assume that the evolution of structure 
between z = 0.1 and z = is negligible. All model galax- 
ies are shifted into redshift space using the distant observer 
approximation. 

Environments are measured using the smoothed top- 
hat number density contrast around each galaxy. As in 
ICroton et all tood ). we define a density defining popula- 
tion (DDP) of galaxies which fix the density contours in 
the redshift space volume containing the model galaxies; 
this is simply all galaxies in the magnitude range — 19 > 
Mbj — 5 logj^p/i > —22. We then consider all galaxies in turn, 
counting the number of DDP neighbours with in a sphere of 
8 /i~^M pc radius, Ng. This scale was found bv lCroton et al.l 
to optimise the simultaneous sampling of both the 
under- and over-dense regions of the survey in a consistent 
way. The mean number density of the DDP is known, p, 
and thus the expected mean number within 8/i~^Mpc is 
also known, Ng = p|7r8^. Hence, the density contrast, Ss, 
for each galaxy is simply 



Pa 



Ng-Ng 
Na 



(2) 



fl=8h-lMpc 

The Millennium Simulation was run using periodic bound- 
ary conditions, and hence we do not need to worry about 
boundary effects in our semi-analytic model galaxy cata- 
logue or Ss calculation. 

To determine the volume fraction that voids occupy in 



the Millennium Simulation we fill the simulation box with a 
large number of randomly placed points, and for each we de- 
termine the number of DDP galaxies within 8 h~^Mpc. This 
provides a uniform sampling of Ss via Equation[2]throughout 
the box. The fraction of random points in a Ss range provides 
the volume filling fraction that this environment occupies. 
This allows us to appropriately normalise the number den- 
sity of galaxies in that environment when calculating each 
luminosity functio n. This m e thod is described in detail in 
Appendix A of iCroton et al.l (|2005l '). 



4 RESULTS AND DISCUSSION 

4.1 The void luminosity function 

We begin with Figure [T] which shows the luminosity func- 
tion of model galaxies for both the population as a whole, 
and all galax ies in void regions. Here we define a void as 
ICroton et all fioOS) did, ^8 < —0.75, i.e. where the density 
within an 8 /i~^Mpc radius of the galaxy is less than 25% the 
mean. Over-plotted ar e the observational r esults for the full 
2dFGRS (filled circles. iNorberg et al.ll2002l') and void g alaxy 
luminosity function (open sq uares , ICroton et al.l l2005h . 

As was the focus of iCroton et a l. (2006), the global 
luminosity distribution of model galaxies is a good match 
to the local observations. Important model ingredients that 
produced this result are the supernova in shallower poten- 
tials that expel disk gas to reduce star formation and flatten 
the faint-end slope, and the inclusion of an AGN heating 
source in large halos to starve massive central galaxies of 
star forming fuel, resulting in the observed bright-end expo- 
nential cut-off, as described in Section (2] Without these two 
critical aspects in the model, the luminosity function instead 
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Figure 3. The Millennium Simulation halo mass function (left panel) and halo mass function in voids (right panel). In both panels 
the halo mass function is broken up into those who host red central galaxies (solid lines) and those that host blue central galaxies 
(dashed lines). Red sequence galaxies occupy the most massive halos in all environments — these halos are subject to the radio mode low 
luminosity AGN heating that ultimately shuts down star formation. Notably, this shutdown begins at approximately the same mass in 
both panels, M^ir ~ 10^^~^^'^ A/q . This implies that the (environment independent) radio mode heating, plus a shift in the halo mass 
function with environment, is sufficient to reproduce the observed abundance of void early-type galaxies seen in Figure [2] 



would have a pow er-law shape reflectin K the underlying mass 
function of halos (|Benson et al.ll2003al ). 

Note that there remain discrepancies in the top curve 
of Figure [2 notably that the model over-predi cts the abun- 
dance of very bright galaxies. As discussed in ICroton et al.l 
l|2006l ). these galaxies represent a population undergoing 
strong star formation and merger-induced starbursts. In 
such ultraluminous infrared galaxies (ULIRGs), nearly all 
the hght from young stars is absorbe d by dust and re- 
radia ted in the mid- to far-infrared (|Sanders fc Mirabell 
1 19961 ) . Improved modelling of the effects of dust are required 
to adequately reproduce the properties of such systems. 

The lower solid line in Figure [T] shows the model result 
for the void galaxy population in the Millennium Simulation. 
The good overall agreement is a result of the physical pre- 
script ions and global parameter choices in the lCroton et al] 
l|2006l ) model and not additional fine-tuning (of which there 
was none) . There is some over-prediction of the very bright- 
est void galaxies. However the discrepancy is not signifi- 
cant enough to expect that improvements to existing as- 
pects of the model, e.g. dust as described above, cannot 
alleviate such differences. Additionally, we find a small over- 
prediction of faint (Mbj — 5 logj^fj/i — 18.5) galaxies. Un- 
fortunately, systematic variations in the observed faint-end 
galax;y lumin osity function ex ist at a level greater t han this 
(e.g. contrast [Cole et al-llioOll to lHuang et al.ll2003h . 

Figure [T] alone answers the challenge posed by Peebles 
and outlined in the Introduction. Specifically, in a ACDM 
Universe the distribution of dark matter halos in voids, cou- 
pled with a realistic model of galaxy formation, is consistent 
with voids as observed in the real Universe. Such regions ar e 
typically not devoid of galaxies as claimed in lPeebleg ([iooj). 

The remainder of this letter will aim to dissect and un- 
derstand the good agreement shown in Figure [l] between 
model and observation. 



4.2 Galaxy colours in voids 

The colour of a galaxy tells us a lot about the relevant 
physics that has been dominant during its evolution. A red 
spectrum typically indicates that a star formation shutdown 
mechanism has been operating for a significant part of the 
galaxy's lifetime. We can use this knowledge, in conjunction 
with our theoretical model, to gain insight into how shut- 
down may occur and to what degree it may (or may not) be 
environment specific. 

To this end, in Figure [5] we again plot the void 
galaxy luminosity function for both the 2dFGRS (symbols, 
ICroton et ai]|2005l ) and semi-analytic model (lines), but now 
broken up by galaxy spectral type (early/late) or colour 
(red/blue). Early- type (triangles) and late- type (squares) 
2dFGRS galaxies are determined u si ng th e principal com- 
ponent analysis of iMadgwick" et al.1 (|2003l ). For the model 
we use the bi-modal galaxy colours to separate red (solid 
line) fro m blue (dashed lin e) at mtj — m^p = 1.07 (see fig- 
ure 9 of ICroton et al.ll2006l ). Comparing model to observa- 
tion. Figure [2] shows agreement (to within 2a) between the 
two sets of early- type/red and late-type/blue void luminos- 
ity functions. This has occurred as a byproduct of matching 
the model to the observed global properties and only these 
- no special void environment physics was needed. Adding 
detail to the physical prescriptions assumed by the model 
would presumably improve the agreement further, but this 
is not our focus here. 

The early-type 2dFGRS void galaxies in Figure[2]return 
our focus to the question asked in the Introduction: how does 
one understand early-type galaxies in voids when void en- 
vironments are typically very gas rich and slowly evolving, 
which tends to promote star formation rather than suppress 
it. At this point it may be valuable to backtrack somewhat 
and revisit the two main mechanisms acting in the model 
that can turn blue galaxies red. The first is important for 
satellite (i.e. non-central) galaxies. Upon infall into a more 
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massive system, any extended hot gas around the (now) 
satelhte is stripped by the denser medium and added to the 
more massive halo. Such "strangulation" drives a satellite 
galaxy to redden rapidly once its remaining cold disk gas is 
exhausted. The second is the radio mode heating discussed 
above and in Section [2] Radio mode AGN operate only in 
central galaxies where the host halo has grown above a crit- 
ical mass threshold, approximately Mvir ~ 10^^'^ Mq. 

So what physical processes have occurred to produce in 
the observed 2dFGRS early-type void population? It may 
be that such red galaxies are simply a population of satel- 
lites and have survived long enough for an effect on their 
colours to be seen. To test this idea using the model we 
perform a simple and revealing exercise. From our knowl- 
edge of which galaxies are central and which are satel- 
lites, we remove those that are satellites and hence those 
that could have experienced the strangulation effect. This is 
shown by the dotted lines in Figure [21 Although satellites 
do contribute to the void early-type luminosity function at 
A'Ihj~5\ogigh ^ — 18.5, they are clearly a minor component 
brighter than this. Indeed, all bright early-type void galaxies 
are centrals in our galaxy formation model. 

The second of our star formation shut-down mecha- 
nisms, i.e. radio mode low luminosity AGN heating, requires 
that central galaxies reside in sufflciently large dark mater 
halos. Such halos are not expected to be common in void 
regions of the Universe where low mass halos dominate. In 
low mass halos the central density of hot gas is not high 
enough to maintain the low Eddington accretion needed to 
power a central AGN outflow. We check this explicitly in 
Figure [S] where we plot the halo mass function for halos in 
all environments (left panel) and void environments (right 
panel) (wide-dashed lines in both). In addition, we break 
the halo mass functions into those halos hosting red central 
galaxies (solid lines) and those hosting blue central galaxies 
(dashed lines). 

Figure [3] reveals the origin of the void early-type pop- 
ulation - they are simply central galaxies that live in halos 
massive enough to have the low luminosity radio mode op- 
erating. In other words, the physics assumed by the model 
that transforms blue galaxies into red operates independent 
of environment but can still produce the observed environ- 
mental trends in voids (note that the transition from blue- 
dominated to red-dominated dark matter halos begins at ap- 
proximately the same mass in both panels, Mvir ~ 10^^ Mq). 
Specifically, it is the shift in halo mass function with envi- 
ronment that changes, not the transformation mechanism 
itself. 

This picture has important consequences for galaxy 
evolution more generally. First, it indicates the necessity 
for active galaxie s in voids, confirmed in a number of 
studies (see, e.g.. IConstantin et"al] |2008| at z ~ 0, and 



iMontero-Dorta et al.l 20081 a,t z ^ 1). Second, it requires 



that void galaxies should have similar properties to those 
in all other environments, assuming the galaxies compared 
are hosted by halos of similar mass. This is a theoretical 
confirmation of numerous observational studies that have 
shown that environments on scales larger than a few Mpc 
(i.e. outside the dark ma tter halo) are not impo rtant for 



galaxy formation (e.g. see lBlanton fc Berlindl 



erences therein). For example, Patiri et al.l 



20071 and ref- 
20061 ) investi- 



gate the properties of void galaxies in the Sloan Digital Sky 



Survey and compare to tho se of a similarly cons tructed mock 
catalogue using the same ICroton et al.l l|2006l ) model used 
here. They found little difference in the mean properties of 
void galaxies relative to the field, specifically in colour, spe- 
cific star formation rate, and morphology. Theoretically one 
may expect large-scale environment to be imp ortant, e.g. 
the e ffect of assembly bias for dark matter halos (|Gao et al.l 
I2OO5!) . However the galaxy population d oes not appear to b e 
overly sensitive to this (see, for example. ICroton et al.ll2007l '). 
A similar conclu sion was recently arr ived at in a complimen- 
tary analysis by [Tinker et~aLl l|200'i1) . 



5 SUMMARY 

It has been suggested that void environments pose a problem 
for galaxy formation theory in a ACDM universe due to the 
apparent over-abundance of dark matter halos in voids rela- 
tive to the observed abundance of galaxies (Peebles 2001). In 
this paper we show that no conflict exists within the current 
observational uncertainty. 

• Our model can reproduce the observed abundance of 
void galaxies, both globally and by colour, without any mod- 
iflcation or additional environment dependent physics. 

• Early-type "red and dead" red sequence galaxies ap- 
pear naturally in the voids. They arises because of a shift 
in the halo mass function in low density environments com- 
bined with an environment independent star formation shut- 
down mechanism efficient above a critical halo mass (here, 
radio mode AGN). Together, these approximately produce 
the correct observed abundance. 

• Some notable consequences follow from our results. For 
example, at a given host halo mass, void galaxies are ex- 
pected to have similar properties on average to those in the 
field, because such galaxies will have had similar evolution- 
ary histories to field and cluster galaxies. 

Voids and void galaxies provide an important probe of 
both cosmology and galaxy formation. Future wide-field sur- 
veys will produce large-scale maps of the Universe out to 
high redshift, from which void evolution can be studied in 
detail. Such work will further constrain galaxy formation 
theory, both in under-dense and the wider field environment. 
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